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Abstract. - For percolating systems, we propose a universal exponent relation connecting the 
leading corrections to scaling of the cluster size distribution with the dynamic corrections to the 
asymptotic transport behaviour at criticality. Our derivation is based on a cluster-resolved scaling 
theory unifying the scaling of both the cluster size distribution and the dynamics of a random 
walker. We corroborate our theoretical approach by extensive simulations for a site percolating 
square lattice and numerically determine both the static and dynamic correction exponents. 
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q , Anomalous, subdiffusivc transport has been widely ob- 
~~ Served in biological systems, e.g., in the cytoplasm [1] and 
I in cell membranes [2] . These findings are attributed to an 
obstructed motion due to the crowded nature of cellular 
■"xj" environments, and they have been related to the subdif- 
~~ < fusivc motion established for percolating systems [3] . The 
^ ^ experiments have yielded a wide spectrum of almost con- 
tinuously changing exponents, but a mechanism gencrat- 
,— ' ing such a continuum of fractal exponents has not been 
identified yet. Although the observation windows cov- 
ered between one and three decades in time, an alternative 
^•interpretation in terms of apparent power laws resulting 
from crossover phenomenona over large windows is con- 
ceivable. 



threshold p c for large clusters, s — > oo. The Fisher ex- 
ponent r can be related to the fractal dimension d{ of the 
infinite cluster via r = l + d/df, provided the space dimen- 
sion d is smaller than the upper critical dimension d c = 6. 
The deviations from the power law for n s are to leading 
order again determined by a power law 



n 3 (p c ) ~ As~ T (1 + Bs- 



for 



(1) 



X 



For a continuum percolation model, it has been shown 
recently that the crossover from pure subdiffusion to nor- 
mal diffusion extends over five decades in time [4] ; in addi- 
tion, the asymptotic behaviour is slowly approached and 
the large corrections cannot simply be ignored. Thus, it 
is of general interest to develop a systematic description 
of universal corrections to scaling in percolating systems. 

In this Letter, we generalise the dynamic scaling theory 
for transport on percolation clusters [5, 6] to include the 
leading corrections. We derive a new exponent relation 
for dynamic corrections that are inherited from the static 
structure. Let us start by introducing some notation and 
a summary of our main results. 

First, the cluster size distribution in the infinite lat- 
tice becomes fractal, n s ~ s~ r , directly at the percolation 



The correction exponent Q > is by renormalisation 
group arguments expected to be again universal, i.e., it 
depends only on the dimension of the system, but not on 
the details of the lattice. Its role is to quantify how fast 
the critical manifold is approached and microscopic de- 
tails become irrelevant. We have determined the value of 
the correction exponent to SI = 0.77 ± 0.04 by extensive 
computer simulations for a two-dimensional (2D) square 
lattice. 

Second, transport of a random walker on the incipient 
infinite cluster at p = p c is anomalous, quantified by the 
power-law dependence of the mean-square displacement 
^ r oo( < ) ~ t 2 !^ for large times t — > oo, where d w denotes 
the walk dimension. The approach to this power law is 
also non-analytic, 

SrKt) ~ A^t 2 ^ (I + Coot-y) for t oo , (2) 

giving rise to the dynamic correction exponent y > 0. 

We shall introduce a cluster-resolved scaling theory for 
a generalised probability distribution incorporating the 
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91/48 a 
4/3 a 
0.77(4) c 



2.530(4) c 
0.875(l) b 
0.64(2) b 



lattice 


continuum 


d w 2.878(l) a 


3.88(3) a 


4.81(2) d 


z 3.036(1) 


5.07(6) 


6.30(3) 


y 0.49(3) c 


0.42(2) c 


0.34(2)° 



Table 1: Static and dynamic exponents for the leading and sub- 
leading critical behaviour. The uncertainty in the last digit is 
indicated in parentheses. Sources: (a) ref. [7], (b) ref. [8], (c) 
calculated from r = 2.186(2) [9], (d) continuum percolation 
theory yields d w = df + 2/u [4, 10], and (e) this work. The 
dynamic universality class does not split for d — 2 [10, 11]. 
The values for z are derived from an exponent relation; the 
predictions for y are based on our proposed relation, eq. (3). 



static and dynamic properties and derive the new expo- 
nent relation 

yd w = fldf . (3) 

By computer simulations, we demonstrate that the as- 
sumptions underlying our scaling hypothesis are indeed 
fulfilled and that the measured correction exponents are 
compatible with the relation above. We have provided ev- 
idence for eq. (3) for d = 3 before [4], and we will test our 
theory here for d = 2; universality allows to choose the 
simple square lattice for this purpose. In Table 1, we have 
compiled some critical exponents for transport on perco- 
lation clusters in two and three dimensions along with the 
predictions from this work. 

Scaling of the cluster size distribution. As the 

probability p for a site to be occupied approaches a certain 
critical value p c , clusters of all sizes emerge and the system 
becomes self-similar. In the infinite lattice, this is reflected 
by the power law tail of the distribution n s ~ s~ r ; in a fi- 
nite lattice of box length L, the largest cluster is expected 
to contain of the order of L d{ sites. Off the percolation 
threshold, the linear extension of the largest finite cluster 
£, also referred to as correlation length, diverges in the 
infinite lattice, £ ~ |e| _ly , where e := (p — p c )/p c denotes 
the reduced distance to the critical point. We now intro- 
duce a scaling assumption for n s that encompasses all the 
singularities as s — > oo, e — > 0, and L — > oo, including also 
an irrelevant scaling field u, 



n s {e,L) = L- d - dl N> 



sL 



-d, 



L 1/v ,uL- 



(4) 



The first scaling variable sL~ d{ states that the size s of 
a cluster should be compared with the size L dt of the 
largest cluster in the box. The second argument sL x l v 
suggests that the system appears similar upon zooming 
in, L i— > XL, provided the correlation length is rescaled 
accordingly, £ t— > A£. The scaling field u is assumed to 
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Fig. 1: Rectification of the cluster size distribution n s with 
the critical power law s~ T yields the amplitude A. Periodic 
(solid, blue) and free (broken, red) boundary conditions are 
compared. The dotted line indicates the leading correction to 
the asymptotic scaling, n s ~ As~ T {l + Bs~ n ). Inset: the de- 
viation from the asymptotic law on a double-logarithmic plot 
follows a power law again with the universal correction expo- 
nent Q = 0.77. 



be the leading irrelevant variable in the renormalisation 
group sense [12], encoding how some microscopic details 
of the system fade out as the critical point is approached. 
Since here the system is finite, n s is an analytic function of 
all its arguments, and since the arguments (s,£, u) appear 
linearly inside the scaling function N, the same property 
holds for N. The leading and next-to- leading behaviour is 
extracted by expanding N to first order in the irrelevant 
scaling field u, 



,{e,L) = L- d - d *UQ{sL- d! ,eL l l v ) 

+ £-rf-*-«Ni(sI,-* eL 1 /'') 



(5) 



with new analytic functions No, Ni, and a factor u has 
been absorbed in A similar form of eq. (5) was already 
given by Margolina et al. [13]. 

Static corrections to scaling. — For large systems, 
L>(, the dependence on the box size may be eliminated, 
and eq. (5) for the cluster size distribution can then be 
written as 

n s {e) = s- T Kt (*r df ) [1 + a'**? (sr*)] (6) 

with new one-parameter scaling functions No (•) for the 
leading behaviour; the superscript ± distinguishes the 
sign of e. The correction is parametrised in terms of 
Nf (•) and the new correction-to-scaling exponent := 
w/df. Specified to the critical point, this yields n s (e = 
0) ~ As~ T (1 + Bs- n ) , with the non-universal amplitudes 
A:=N±(0) and B := 3\ff (0). 

A measurement of the correction exponent Q requires 
high quality data for the cluster size distribution n s . The 
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deviation from leading order has to show an extremely 
small signal-to-noise ratio in order to follow the power law 
sufficiently long. This implies a precise knowledge about 
the leading order, n s ~ As~ T , including the prefactor A. 
We have simulated 169,000 square lattices at the thresh- 
old, 1 p c = 0.592746 [14], using the multi-labelling algo- 
rithm by Hoshcn and Kopclman [15]. We have used peri- 
odic boundary conditions and large lattices with a linear 
extend of L = 45,000 sites to minimise the effects due to a 
finite simulation box [16]. Figure 1 shows a rectification of 
the obtained distribution n s : the critical power law is fol- 
lowed over more than five decades in s. At small cluster 
sizes, the critical law is approached with another power 
law, s~ n , as expected. The figure also displays results for 
free boundary conditions, where pronounced finite-size ef- 
fects, however, render a precise observation of the leading 
critical behaviour nearly impossible. 

In two dimensions, we take advantage of the Fisher 
exponent being known exactly, r = 187/91. We have 
determined the amplitude of the leading order to A = 
0.016921(4) by linear regression over various intervals in 
10 4 < s < 10 7 . Over these three decades, our data deviate 
from the leading power law by less than l%o. The devia- 
tion from the asymptotic law is shown in the inset of fig. 1 
for smaller cluster sizes, s < 10 5 . On the double logarith- 
mic plot, the data nicely follow a straight line again over 
more than two decades in the cluster size. But even with 
such accurate data, the precise extraction of the exponent 
is quite intricate. Since we need an asymptotic fit for large 
(but not too large) s, we varied the fit interval for the lin- 
ear regression in 10 2 < s < 10 4 ; this procedure, however, 
seems to underestimate the exponent slightly in our case. 
A more reliable approach for an asymptotic fit is to do 
another rectification, now of the correction: the best con- 
vergence of s Q (n s /As~ T — 1) — ► B was found for ft = 0.77 
with the above value of A, yielding B = —1.01. The result 
slightly depends on A; varying the amplitude within the 
above uncertainty, we have obtained error bounds on the 
correction exponent, = 0.77 ± 0.04. 

Universality of fl suggests it is independent of the lat- 
tice. Indeed, our result is fairly compatible with the result 
for the triangular lattice [17], where was found to be be- 
tween 0.71 and 0.74; we attribute the discrepancy to the 
limited observation window for the correction, s < 10 3 , 
22 years ago. 

Dynamics of a random walker. Transport on 
percolating systems is expected to be anomalous since a 
random walker ( "the ant" ) has to explore the self-similar 
structures of the ramified clusters. The frozen disorder 
leads to subdiffusive motion on the infinite cluster directly 
at the percolation threshold. An independent critical ex- 
ponent d w > 2, referred to as walk dimension, charac- 
terises the restricted mean-square displacement (MSD) de- 

lr The simulations required a total computation time of 90 hours 
on 100 cores of the AMD Opteron 285 2.6 GHz processor with 2 GB 
of RAM per core, which limited the size of the lattices. 



fined as Sr^t) := ([R(t) - R(0)} 2 \ R(0) £ 6^) ~ t 2 /^, 
where the average is taken over different realisations of 
the disorder and R(t) denotes the trajectory of the ant; 
some values of d w are listed in Table 1. If walkers 
on all clusters arc considered, the exponent is modified 
since the average of clusters of all sizes is weighted with 
the fractal cluster size distribution [18]; the unrestricted 
MSD follows <5r 2 v (i) ~ t 2 l z with the dynamic exponent 
z := 2d w /(2 + d[ — d) > d w . Away from the critical point, 
the long-time behaviour is either diffusive, 5r 2 v (t) ~ 2dDt, 
for p > p c or localised, Sr 2 v (t — > oo) = £ 2 , below the 
transition, p < p c . The diffusion coefficient D and the 
localisation length £ exhibit power law behaviour in e. 

Cluster- resolved scaling. The statistical information 
of the transport dynamics is encoded in the van Hove self- 
correlation function G(r,t) :~ ((5(i?(t) — R(0) — r)), i.e., 
the probability that the ant has been displaced by a vec- 
tor r in time t. Here we suggest a cluster-resolved scaling 
theory for the van Hove function including the leading cor- 
rection to scaling; in particular, we will derive a relation 
between static and dynamic exponents. The principle idea 
is that at the critical point all clusters resemble each other 
up to proper rescaling and therefore induce similar dynam- 
ics after readjusting the time scale. To take advantage of 
the self-similarity, we generalise the van Hove function to 
also resolve the cluster size, 

G s (r,t):=(8(R(t)-R(0)-r)\R(0)ee s ) , (7) 

where C s denotes a cluster of size s < oo. After averaging, 
the information on the detailed structure of the cluster is 
lost, and we keep only its size as indicator for the dynam- 
ics. The ant is trapped for finite s, whereas for s = oo 
the infinite cluster allows for long-range transport. Let us 
introduce also the joint probability P s that the ant is on 
a cluster of size s and has moved a vector r in a given 
time i, as well as a corresponding probability for the 
infinite cluster [19], 

P s (r,t;e) := sn s (s) G s (r,t;e) , (8a) 
P 00 (r,t;e):=P(e)G 00 (r,t;e); (8b) 

P(e) is the probability that a given occupied site belongs 
to the infinite cluster. We now postulate a scaling relation 
for the joint probabilities in the same spirit as for eq. (4), 

P s (r,t;e,L) = 

L- 2d P F (r/L,iiT^ S L-* eL^uir") , (9a) 
Poc(r, t; e, L) = 

L d t -2d p oo (j/L,tL- d ™ ,eL l f v ,uL- u ) , (9b) 

where the finite size L of the system has been explicitly 
considered. (For a finite system, the infinite cluster is de- 
fined as the largest cluster of a given percolating lattice.) 
We have assumed that the leading irrelevant scaling vari- 
able u is the same as in the static case. Since the system 
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is finite, Pp and are analytic functions in all their ar- 
guments. The choice of the arguments is motivated by 
the well-known scaling relations for the linear extent of 
walks (exponent d w ), the cluster size (<if), and the corre- 
lation length (y). The prcfactor ensures consistency with 
the scaling form for the cluster size distribution, eq. (4). 
Since the van Hove functions G s arc normalised, sn s is ob- 
tained by integrating eq. (8a) over all r, compatible with 
the scaling hypothesis for n s . Integration of eq. (9b) yields 
the familiar result for the strength of the infinite cluster, 
P ~ le^, with exponent ft = v(d — df) as L — ► oo. To sim- 
plify notation, we suppress the dependence of the scaling 
functions on u in the remainder of this section. 

The limit of infinite system size eliminates the depen- 
dence on L, 

P s (r,t;e) = R- 2dt 5 , ¥ {r/R s ,tR- d ™,eR 1 J v ) , (10a) 
Poo(r,t;e) - e*" 2 ^ (r/U^) , (10b) 

where in the new scaling function CPp, we have replaced the 
cluster size s by the linear extent of s-clusters, R s ~ s 1 ^ 1 , 
which serves as reference length scale. Note that CPp is still 
analytic in all arguments. For the van Hove functions, this 
implies a scaling form 

G s (r,t;e) = R- d S F (r/ R s ,tR; d * ',eR\'") , (11a) 
Goo{r, t; e) = C d 5oo (r/£, tr d ~) ■ (Hb) 

The second moments define the cluster-resolved MSDs and 
possess the scaling forms 

5r 2 s (t;e) = ^ d *S^(tRj d ',eRy v ) , (12a) 
Srl(t 1 e) = t 2 / d ~SKl(tC d ~) ■ (12b) 

Note, that most clusters are much smaller than the cor- 
relation length, R s <C £, implying |e|i?s^ <C 1. For these 
clusters, Sr 2 (t; e) is insensitive to the distance to the crit- 
ical point e, and thus eq. (12a) essentially constitutes a 
one-parameter scaling prediction, which we will discuss 
below. The second relation, eq. (12b), is the familiar re- 
sult for a random walker restricted to the infinite cluster. 

We have generated a total of more than 1000 trajectories 
on 20 square lattices with a linear extent of L = 8,000 sites 
for several densities p each. Periodic boundaries were used 
to minimise finite-size effects, the starting points of the ant 
were chosen at random, and the size of the clusters was 
determined using a union-find algorithm [20]. 

The cluster- resolved MSDs Sr 2 (t) at the threshold are 
exhibited in fig. 2; all of them follow the same power 
law f 2 / dw before they saturate. The scaling prediction 
eq. (12a) implies that the respective long-time limits yield 
the radii of gyration and grow as R s ~ s 1 ' dl . The crossover 
time from subdiffusion to localisation is expected to scale 
as t x ~ i?f w . Figure 3 demonstrates that our data exhibit 
this scaling remarkably well; they bridge the crossover re- 
gion from anomalous transport to localisation, which ex- 
tends over three orders of magnitude in time. We have also 
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Fig. 2: Cluster-resolved MSDs 5r 2 a {t) at the threshold (solid 
lines); Sr^if) is shown for several densities above the threshold. 
Distances are measured in units of the lattice spacing and time 
refers to the number of hopping attempts. 

included the MSD on the infinite cluster Sr 2 ^ (t) above the 
threshold in fig. 2. The crossover to diffusive transport 
occurs once the MSD is of the order of £ 2 , i.e., right at 
the length scale above which the system looks homoge- 
neous, or in terms of times, t x ~ £ dw , in consistency with 
eq. (12b). The scaling prediction is illustrated in fig. 4, 
corroborating an earlier study (fig. 6.2 in ref. [5]). In sum- 
mary, the observed scaling behaviour validates the notion 
of self-similarity also for the structures of the clusters and 
the transport on them, provided the system is close to the 
critical point. 

All-cluster-averaged MSD. The all-cluster-averagcd 
MSD 5r 2 v (t) is calculated from the data for the cluster- 
resolved MSDs and the weights sn s and P, 

Srl v (t;e) =J2^s(e)Sr 2 (t) + P(s)Srl(t) . (13) 

s 

Using the scaling predictions eqs. (6) and (12), one easily 
derives Sr 2 v (t;e) = i 2 / z <53? 2 t (t^~ dw ), where the dynamic 
exponent z emerges. At criticality, the infinite cluster 
has zero weight and the finite clusters contribute an in- 
finite hierarchy of cluster sizes, resulting in an unbounded 
MSD, diverging as Sr 2 v (t) ~ t 2 l z . Below the threshold, 
Sr 2 (t — » oo) ~ R 2 implies that the long-time limit of 
the averaged MSD measures the mean-square cluster size, 
^ r av(^ — * °°) =: where £ ~ \e\~ u+l3 / 2 . Consistency 
with the dynamic scaling for 5r 2 v (t) requires £ z ~ £ dw , 
allowing to write the scaling form more transparently, 

5r 2 av (t) = t 2 ^6X 2 ± (U- z ) . (14) 

On large time scales when 5r 2 v (t) ^> £ 2 , merely the infinite 
cluster contributes to the MSD and transport becomes 
diffusive. With the ansatz 5r 2 v (t) ~ Dt, eq. (14) dictates 
that the diffusion coefficient D vanishes at the threshold 
as D ~ \e\^ where [i — v(d w — 2) + ft. 
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Fig. 3: Cluster-resolved data from fig. 2 rescaled with s 1 ' df , 
which is proportional to the radius of gyration R a . 



Close to the percolation threshold, our numerical results 
for the MSD 5r 2 v {t\e) follow the critical law, 5r% v (t) ~ 
t 2 / z , until the curves fan out to either diffusive or localised 
behaviour (inset of fig. 4). We have rescaled the simulated 
data according to eq. (14) and find excellent data collapse, 
see fig. 3. The crossover region covers approximately four 
nontrivial decades in rescaled time. Let us emphasise that 
to achieve a data collapse of similar quality for the 3D 
Lorentz model, it was necessary to consider explicitly the 
leading corrections to scaling [4]. 

Corrections to dynamic scaling. The leading dy- 
namic corrections to scaling are derived from the scaling 
form of the joint probabilities for finite and infinite clus- 
ters, eqs. (9), along the same lines as for the cluster size 
distribution. For the MSD on the infinite cluster this im- 
plies 5r£,(i;e) = t 2 / d »(5R^(^- d -,<- w ). Expanding to 
first order in the irrelevant scaling variable u yields 

Srl(t) = t^SXl {tC d ~) [1 + t-VA^ (tC d -)] , (15) 

where the correction is quantified by the new exponent 
y = uj/d w and a new scaling function A OQ (-). Specifically 
at the critical point, the cluster-resolved scaling hypothe- 
sis predicts Sr^t) = A^t 2 / 4 ™ (1 + C^y) with the non- 
universal amplitudes A^ := <53^(0) and Coo := A M (0), 
compare eq. (2). Thus, we have derived the exponent re- 
lation, yd w = cu = fldf, between the dynamic and static 
correction exponents y and f2, mentioned in the introduc- 
tion. 

A similar chain of arguments for the MSD on the finite 
clusters extends eq. (12a) at criticality to 

Sr 2 s (t) = t 2/d -SJl 2 F (tR- d -) [1 + t- y A F (tR; d ")] . (16) 

Observing that the MSDs on a finite and the infinite clus- 
ter are indistinguishable at not too long time scales, t -C 
Ef-, fS we infer £Rf (0) = ^oo and A F (0) = C^. Calcu- 
lating the all-cluster average, eq. (13), from eqs. (6), (15), 



Q 



10" 



10' 



10 u 



10" 



10"' 





























5R 2 + (i) ' diffusion 


subdiffusion 














10 8 


5r 2 (t) 






10 6 


\ \ localisation 


" 10 4 






\ i~ 2/z 


10 2 




Time t 






10 2 10 4 


10 6 10 8 1C 


10 \ 



10' 



10 u 



e 
o 

0.001 

0.003 

0.005 

0.01 

0.02 

0.03 

0.05 

0.08 

-0.001 

-0.003 

-0.005 

-0.01 

-0.02 

-0.03 

-0.05 

-0.08 



10"' 



10"' 



10"' 



10"' 



10 u 



10" 



1(T 



Rescaled time t = t\e\ 



tt 



Fig. 4: All-cluster-averaged MSDs 5r 2 v (t; e) close to the thresh- 
old. While the inset shows the raw data, rescaling according 
to eq. (14) collapses all data on the scaling functions £R±Q. 
Right axis: rescaled MSDs on the infinite cluster, 5r^o(t;e), 
from fig. 2, following eq. (12b). 



(16) leads to 
Sr 2 av (t) = 



t 2 ' z 8R 2 ± 



'■) [i + t-»A±(tr*)] , (17) 



where the scaling variable can equivalently be written as 
££ -dw oc tt~ z . We have suggested such a form already for 
the 3D Lorentz model [4], and here we provide the basis 
for this ansatz. In particular, we shall test the prediction 
at the critical point, 



5rl v (t)~A av t 2 ' z (l + C av t-y) for t 



(18) 



For a numerical analysis of the dynamic corrections to 
scaling, we have simulated 4,070 trajectories of 10 11 steps 
on periodic square lattices with a linear extent of L = 
90,000 sites at the critical density; 2 the starting points 
were chosen at random on arbitrary clusters. To increase 
computational efficiency, we utilised the statistical inde- 
pendence on the bit-level of the random number gener- 
ator GFSR4 [21]: calculating one random 32-bit integer 
allowed the realisation of 16 hopping attempts. Figure 5 
shows the rectified MSD 5r 2 v {t)/A m t 2 / z at the threshold; 
a fit for the amplitude yields A^ = 1.204(1). The leading 
asymptotic behaviour is very well captured over 6 decades 
in time with the exponent z obtained from conductivity 
measurements, see table 1. On a double- logarithmic scale, 
the deviations from the leading behaviour are compatible 
with a power law decay with exponent y = 0.49. There is 
still a relatively large, overall statistical error in our data 
of about 1% since the heterogeneous average over clusters 
of all sizes converges quite slowly. Better convergence is 



2 Each trajectory required a computing time of about 19 hours on 
a single core of the AMD Opteron 8218 2.6 GHz processor. 
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Finally, our analysis suggests to interpret experiments 
on anomalous diffusion in disordered environments in 
terms of universal exponents, including leading corrections 
to the asymptotic behaviour, rather than fitting with phe- 
nomcnological exponents. Percolation provides one possi- 
ble universality class for such an approach. 
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Fig. 5: Rectified MSDs restricted to the infinite cluster (solid, 
blue) and averaged over all clusters (broken, red) at the critical 
density. Data for 5rj^,(t) are rescaled with Acct 2 ^ d ™ and for 
Sr 2 v (t) with A av t 2 ^ z . Inset: a double-logarithmic plot of the 
relative deviations from the asymptotic behaviour shows that 
the leading corrections indeed follow the predicted power law. 
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derivation relies on a cluster-resolved scaling theory unify- 
ing structure and dynamics, which has been corroborated 
in detail by simulations for the square lattice. 

The corrections to scaling discussed here have their ori- 
gin in the same leading irrelevant scaling field u for the 
statics and dynamics. Yet, the motion of a random walker 
on the fractal clusters introduces new, genuinely dynamic 
scaling fields, one relevant and infinitely many irrelevant 
ones. In our approach, we have anticipated that the lead- 
ing irrelevant scaling field is still provided by the static 
field u, in agreement with our simulations. Analytic de- 
formations of the scaling fields typically lead to analytic 
corrections with integral relative exponents [12]. Here, we 
find £1, y < 1 and conclude that the universal corrections 
dominate over the analytic contributions. 
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